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ABSTRACT 

Because  underwater  gravity  surveys  are  usually  conducted  in  coastal  waters, 
several  unique  corrections  are  required  to  reduce  the  observed  data  for  calculation 
of  free-air  and  Bouguer  anomalies.    The  most  logical  sequence  for  the  application 
of  corrections  would  be:    terrain  correction  (if  necessary)  to  reduce  bathymetric 
irreqularities  to  a  flat,  horizontal  bottom;    removal  of  the  upward  attraction  of 
the  water  above  the  meter;    free-air  correction  to  bring  the  meter  measurement  to 
the  surface  of  the  reference  spheroid;    filling  in  the  ocean  below  the  spheroid  with 
rock  for  a  Bouguer  anomaly  calculation  or  with  water  for  a  free-air  anomaly 
calculation.    The  terrain  correction  must  be  handled  carefully  as  the  bathymetric 
irregularities  are  immersed  in  water  of  finite  density  and,  due  to  the  close 
proximity  of  land,  topographic  irregularities  immersed  in  air  may  be  within  the 
range  of  the  correction  graticule  grid.    Graticule  compartment  weighting  factors 
for  making  this  correction  may  be  readily  determined  by  breaking  the  bottom- 
land surface  profile  into  several  segments  delimited  by  the  c,  avimeter  depth  and 
the  sea  surface . 
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CORRECTIONS  FOR  UNDERWATER  GRAVIMETRY 


INTRODUCTION 

The  concept  that  present  shorelines  represent  generally  insignificant 
discontinuities  to  the  stratigraphy  and  structure  of  crustal  rocks  found  on 
land  and  the  adjacent  continental  terrace  is  basic  to  geology.   However,  the 
problems  of  safe  ship  navigation,  sampling  outcrops  and  crustal  rocks  buried 
by  sediments,  and  conducting  continuous  geophysical  surveys  in  the  shallow 
coastal  waters  present  difficulties  in  the  testing  of  this  continuity.   For 
example,  the  correlation  of  land  gravity  measurements  with  sea-surface  gravity 
surveys  often  requires  interpolation  of  data  across  5  to  10  km  of  unsurveyed 
terrain  between  the  shoreline  and  shiptracks.   One  technique  for  filling  this 
void  is  through  the  use  of  the  underwater  (bottom)  gravimeter  from  a  shallow- 
draft  vessel.   Although  underwater  gravimetry  requires  laborious  field  work, 
the  fact  that  the  meter  is  placed  on  the  ocean  bottom  increases  the  accuracy 
of  the  observed  values  over  those  measured  at  the  sea  surface  by  avoiding  the 
problems  of  removing  the  ship  motion.   Also,  in  coastal  waters  sea-surface 
gravity  data  may  be  unusable  due  to  the  degredation  caused  by  severe  ship  mo- 
tion and  excessive  course  changes  to  avoid  navigational  hazards. 

Although  most  discussions  of  underwater  gravity  surveys  include  some  men- 
tion of  data  corrections,  there  is  a  general  dearth  of  detail  for  a  person  con- 
fronted by  this  situation  for  the  first  time.   In  addition,  the  terrain  correc- 
tion for  bottom  gravity  data  can  become  an  important  factor  along  continental 
shelves  cut  by  steep  submarine  canyons  and  valleys;  the  application  of  this 
correction  to  data  collected  on  the  ocean  bottom  raises  some  important  pro- 
blems to  the  investigator.   The  object  of  this  paper  is  to  present  a  simpli- 
fied discussion  of  data  reduction  applicable  to  underwater  surveys.   Tradi- 
tionally, gravity  corrections  are  discussed  in  the  following  order:   free-air, 
Bouguer,  and  terrain.   This  order  is  significant  in  terms  of  increasing  diffi- 
culty of  calculation,  but  it  is  not  necessarily  the  most  logical  order  for 
understanding  the  physical  significance  of  these  corrections  and  their  proper 
application  to  a  specific  geologic  region. 

Although  underwater  gravity  surveys  have  been  conducted  from  a  diving 
bell,  lowered  to  the  ocean  bottom  with  a  gravimeter  and  its  operator  inside 
(Frowe,  1947),  the  most  practical  method  employs  a  gravimeter  encased  in  a 
pressure  housing  and  controlled  remotely  from  the  ocean  surface  through  elec- 
trical conductors  in  the  lowering  cable.   Pepper  (1941)  and  LaCoste  (1967) 
discussed  the  functioning  of  such  gravimeters  and  Beyer  et  al.  (1966)  deter- 
mined the  feasibility  and  precision  of  measurements  made  with  a  modified 
LaCoste  and  Romberg  gravimeter  in  water  depths  up  to  about  900  m.   At  the 
Naval  Postgraduate  School,  the  LaCoste  and  Romberg  Model  H6G  Underwater  Gravi- 
meter is  used  for  shallow-water  surveys.   This  instrument  is  essentially  a 
land  meter  mounted  inside  a  shell  of  two  thick  aluminum  hemispheres  supported 


by  a  triangular  base.   Regulation  of  the  temperature  of  the  meter  is  accom- 
plished with  a  thermister-transistor  circuit  in  the  meter  housing  within  the 
pressure  shell.   Clamping  and  unclamping  of  the  gravimeter  beam  or  weight, 
leveling  of  the  meter  within  the  shell,  and  manipulation  of  the  measuring 
screw  are  all  remotely  controlled  from  the  sea  surface.   The  maximum  bottom 
slope  permissible  for  remote  leveling  of  the  meter  is  15°.   From  the  ship,  the 
meter  is  lowered  and  raised  by  a  hydraulic  winch  containing  a  slip  ring  assembly 
for  electrical  continuity  between  the  meter  and  the  control  box.   The  control 
box  contains  functions  for  meter  temperature  control,  leveling,  beam  clamping 
and  unclamping,  and  measuring  screw  control.   Also  on  the  control  box  are  indi- 
cators showing  possible  flooding  of  water  into  the  meter,  unacceptable  tilting 
of  the  meter  due  to  excessive  bottom  slope  or  tipping  of  the  meter  on  the  bottom, 
and  the  water  depth  as  measured  by  a  pressure  sensor  on  the  meter,  calibrated 
for  salt  water.   In  50  m  of  water,  the  time  required  for  a  measurement  is  approxi- 
mately 4  minutes.   The  accuracy  of  the  meter,  determined  by  periodic  checks  at 
a  base  station,  is  between  0.1  and  0.2  mgal  (La  Coste,  1967;  Brooks,  1973),  com- 
pared with  accuracies  of  1  to  2  mgal  possible  for  sea  surface  gravimeters  when 
good  navigational  control  is  available  and  a  reasonably  calm  sea  state  prevails 
(Wall,  Talwani,  and  Worzel,  1966;  Beyer  et  al. ,  1966:  Dehlinger  and  Chiburis, 
1972). 

CORRECTION  MODELS 

To  demonstrate  the  various  corrections  to  be  applied  to  the  observed  under- 
water gravity  values,  the  following  models  will  be  described:   flat,  horizontal 
bottom:  smooth  sloping  bottom;  irregular  bottom.   General  equations  for  the 
corrections  are  given  along  with  specific  formulas  assuming  densities  for  sea 
water  and  crustal  rocks.   These  formulas  may  be  modified  for  particular  loca- 
tions using  more  appropriate  densities. 

Flat,  horizontal  bottom 

The  simplest  model  to  begin  with  in  discussing  underwater  gravity  correc- 
tions is  that  of  a  flat,  horizontal  bottom  at  depth  h  below  the  sea  surface 
and  a  difference  between  the  sea  surface  and  the  reference  spheroid  (approxi- 
mated by  mean  sea  level)  of  a  vertical  distance  (h  -  e)  (a  tidal  correction) 
(Fig.  1).   In  this  model  no  terrain  correction  is  necessary  due  to  the  absence 
of  bathymetric  relief. 

The  first  step  in  correcting  the  observed  gravity  is  to  remove  the  gravi- 
tational effect  of  the  water  above  the  meter.   This  is  accomplished  by  calcu- 
lating the  effect  of  a  'Bouguer  plate'  of  water  of  density  a     and  thickness  h. 
The  general  equation  for  this  correction  (Ag-)  is: 

Ag.  =  2ira  Gh,  (D 

°1      w 

—8 
where  G  is  the  universal  gravitational  constant  (6.670  x  10   cgs  units). 
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For  a  water  density  of  1.02  7  g/cm  and  h  measured  in  meters,  eq.  1  reduces  to: 

Ag  =  (0.0430  h)  mgal.  (2) 

Since  the  water  is  attracting  the  gravimeter  beam  or  weight  upward,  this  cor- 
rection is  added  to  the  observed  gravity. 

The  second  step  in  the  correction  requires,  in  essence,  bringing  the  meter 
from  the  ocean  bottom  to  the  reference  spheroid  using  the  free-air  correction. 
The  general  equation  for  this  correction  (Ag  )  is: 

Ag2  =  2  GMe/R3,  (3) 

where  e  is  the  depth  of  the  bottom  below  the  reference  spheroid  (mean  sea 
level)  (Fig.  1),  Mis  the  mass  of  the  Earth  (5.976  x  1027  g) ,  and  R  is  the 
Earth's  radius  (6371  km)  (MacDonald,  1966),  or: 

Ag2  =  (0.3083  e)  mgal,  (A) 

where  e  is  measured  in  meters.   This  correction  is  subtracted  from  the  observed 
gravity  as  the  meter  is  always  below  the  reference  spheroid.   The  value  for 
(h  -  e)  is  numerically  positive  for  tidal  heights  above  mean  sea  level  and 
negative  for  heights  below  this  level. 

The  third  and  final  step  requires  filling  in  the  ocean  (now  air  between 
the  reference  spheroid  and  the  bottom)  with  rock  if  these  gravity  observations 
are  to  be  correlated  with  land  measurements.   This  step  is  numerically  analo- 
gous to  the  'Bouguer'  correction  for  sea-surface  gravity  surveys.   The  equa- 
tion for  this  correction  (Ag~)  is : 

Ag3  =  27TarGe,  (5) 

where  a    ,  the  rock  density  for  the  Bouguer  plate,  is  often  taken  to  be  2.67 
g/cm^  (Dobrin,  1960).   Grant  and  West  (1965,  p.  240)  provide  a  discussion  on 
the  proper  selection  of  the  value  for  the  Bouguer  plate  density  as  it  is  impor- 
tant to  relate  this  value  to  the  local  geology  in  order  to  obtain  a  good  corre- 
lation with  land  gravity.   Using  a  rock  density  of  2.67  g/cm^  (a  value  that  will 
be  used  throughout  this  paper),  eq.  5  reduces  to: 

Ag3  =  (0.1119  e)  mgal,  (6) 

where  e  is  measured  in  meters.   This  correction  is  added  to  the  observed  gra- 
vity as  mass  is  being  added  below  the  reference  spheroid. 

Combining  all  these  corrections  into  one  expression  (Ag  )  gives: 


Agb  -  Ag1  -  Ag2  +  Ag3 

=  2Tra  Gh  +  2G(tto"  -  M/R3)e,  (i) 

w         r  v ' ' 


or 


Ag  =  (0.0430h  -  0.1964e)  mgal.  (8) 

Since  the  tidal  height  (h  -  e)  will  be  generally  small  relative  to  h,  this  total 
correction  will  usually  be  negative. 

Because  of  the  absence  of  bottom  relief  for  this  model,  eliminating  the 
necessity  for  a  terrain  correction,  the  Bouguer  anomaly  (BA)  is  defined  by  the 
relationship : 

BA  =  gQ  +  Agb  -  gt,  (9) 

where  g   is  the  gravity  value  observed  by  the  underwater  meter  and  corrected 
for  earth  tides,  instrument  drift,  and  earth  curvature,  and  tied  to  a  reference 
point  within  the  international  network  of  absolute  gravity  values.   The  term  g  , 
the  theoretical  value  for  gravity  on  the  surface  of  the  reference  spheroid,  is 
given  as  a  function  of  latitude  (X)  by  the  International  Gravity  Formula 
(Heiskanen  and  Vening  Meinesz ,  1958): 

g  =  [978.0490  (1  +  0.0052884  sin2X  -  0.0000059  sin22X)]  gal.     (10) 

Bouguer  anomalies  as  calculated  from  eq.  9  for  underwater  gravity  stations  may 
be  correlated  with  Bouguer  anomalies  for  land  stations. 

Sea-surface  gravity  measurements,  when  corrected  for  horizontal  and  vertical 
ship  accelerations,  cross-coupling  effects  due  to  interaction  between  horizontal 
and  vertical  accelerations  on  the  meter,  instrument  drift,  earth  and  ocean  tides, 
earth  curvature,  and  the  Eotvos  effect,  give  free-air  gravity  anomalies  when  com- 
pared with  g   as  formulated  in  eq.  10.   In  order  to  correlate  these  anomalies 
with  underwater  gravity  measurements,  a  slightly  different  set  of  corrections  to 
the  observed  bottom  value  is  necessary.   The  first  step  is  to  remove  the  upward 
attraction  of  the  overlying  water  using  Ag   from  eq.  1.   Next,  the  meter  reading 
is  reduced  from  the  bottom  value  to  the  value  on  the  reference  spheroid  through 
use  of  Ag   from  eq.  3.   Finally,  the  ocean  is  again  filled  with  water  below  the 
reference  spheroid  using  the  following  correction  (Ag, ) : 

Ag.  =  2-n-a  Ge ,  (11) 

4      w 


or: 


Ag4  =  (0.0430  e)  mgal,  (12) 


where  e  is  measured  in  meters.   This  correction  is  positive  as  the  mass  of 
water  added  increases  the  downward  gravitational  attraction  on  the  meter  beam 
when  located  on  the  reference  spheroid.   The  combination  of  Ag  and  Ag   is  called 
the  double  Bouguer  correction  (Sheriff,  1973,  p.  61).   Combining  these^correc- 
tions  into  one  expression  (Ag  )  gives  the  relationship: 


'f 


Agf   =   Ag1   -   Ag2   +   Ag4 


=   2ttct   Gh  +  2G(tto      -  M/R3)e.  (13) 

w  w  v-l-j; 

Using  values  for  a    ,  G,  M,  and  R as  defined  previously  and  measuring  h  and  e 
in  meters,  this  equation  reduces  to: 

Agf  =  (0.0430h  -  0.2652e)  mgal.  (14) 

As  with  the  expression  for  Ag   in  eq.  8,  this  correction  will  usually  be  negative. 
Thus,  the  free-air  anomaly  (FAA)  for  the  underwater  gravity  measurement  correc- 
table with  sea-surface  free-air  anomalies  is  defined  as: 

FAA  =  gQ  +  Agf  -  gt.  (15) 

Smooth  sloping  bottom 

In  this  model  (Fig.  2),  the  terrain  correction  is  introduced  to  reduce  the 
underwater  gravity  measurements.   Although  the  bottom  is  sloping,  there  are  no 
bathymetric  or  topographic  irregularities  in  the  model.   The  most  logical  approach 
to  this  problem  is  first  to  make  terrain  corrections  that  result  in  a  model  for 
a  flat,  horizontal  bottom,  and  then  to  reduce  the  gravity  data  in  the  same  manner 
as  was  discussed  in  the  previous  section.   Because  the  bottom  is  immersed  in 
water  of  finite  density  rather  than  in  air  with  a  negligible  density  and  because 
some  of  the  'bottom'  may  be  above  the  reference  spheroid  due  to  the  close  proximity 
of  the  shoreline,  certain  modifications  of  the  weighting  factors  for  the  zones 
in  the  terrain  correction  graticule  are  necessary.   The  terrain  correction  chart 
and  tables  developed  by  Hammer  (1939)  may  be  used*  these  tables  assume  a  density 
of  2.0  g/cm^  for  the  material  making  up  the  relief.   Use  of  digital  computers 
in  place  of  the  charts  and  tables  can  significantly  reduce  the  time  required  to 
make  these  corrections  (Kane,  1962;  Bott,  1959). 

The  first  step  in  applying  the  terrain  correction  to  the  underwater  station 

in  Fig.  2  is  to  fill  in  segment  A  (a  wedge  of  water  between  the  station  depth 

and  the  bottom  seaward  of  the  station)  with  rock  of  density  a    .   Since  water  of 

density  a   already  fills  this  wedge,  the  correction  density  or  the  filled  wedge 

has  the  value  (a   -  a  ).   Assuming  densities  of  1.027  g/cmj  and  2.67  g/cm^  for 

a     and  a    ,  respectively,  and  using  Hammer's  correction  tables,  the  total  terrain 

correction  for  compartments  in  this  segment  is  multiplied  by  a  factor  of  0.822 

[i.e.,  (a      -a   )/a    ,  where  a      is  the  density  used  to  make  up  the  selected 
'   r    w   t        t 


correction  table] .   This  correction  is  added  to  the  observed  gravity  as  it  adds 
mass  below  the  meter  depth. 

The  next  step  requires  removing  the  effects  of  the  wedge  of  rock  above  the 
station  depth  (segments  B  and  C  in  Fig.  2).   There  are  two  parts  to  this  step 
because  the  material  in  segment  B  is  immersed  in  water,  while  in  segment  C  it  is 
immersed  in  air.   The  total  terrain  correction  for  segments  B  and  C  is  first 
calculated  assuming  a  density  of  a    ;  again  using  Hammer's  tables  and  assuming  a 
value  of  2.67  g/cm  for  o    ,  this  correction  must  be  multiplied  by  a  factor  of 
1.335  (i.e.,  a   /a   ) .   This  part  of  the  correction  is  also  added  to  the  observed 
gravity  as  it  represents  removal  of  mass  above  the  meter  location.   At  this  point 
it  is  obvious  that  too  much  mass  has  been  removed  by  this  correction  because  in 
segment  B  part  of  the  terrain  effect  of  the  rock  is  balanced  by  water  of  finite 
density.   Thus,  the  second  part  of  this  step  of  the  terrain  correction  involves 
calculating  the  effect  of  filling  in  segment  B  with  water.   The  total  terrain 
correction  for  this  segment  must  be  multiplied  by  a  factor  of  0.514  (i.e., 
a  /a  )  if  using  Hammer's  tables.   This  part  of  the  correction  is  subtracted  from 
the  observed  gravity  as  it  represents  the  addition  of  a  slab  of  water  of  thickness 
h  (segment  B)  above  the  meter  depth. 

In  summary,  the  terrain  correction  involves  filling  in  segment  A  with  rock 
rather  than  water,  filling  in  segment  B  with  water  rather  than  rock,  and  re- 
moving the  effect  of  the  rock  above  the  sea  surface  in  segment  C.   For  the  simple 
model  demonstrated  by  this  case,  Sandberg  (1958)  gave  tables  of  terrain  correc- 
tions for  an  inclined  plane  with  slope  angles  of  1°  to  30°  based  on  Hammer's 
graticule  zones  and  density.   If  the  actual  bathymetry  and  topography  of  the 
study  area  can  be  approximated  by  a  uniform  slope,  these  tables  may  be  used  as 
outlined  above  for  Hammer's  tables. 

The  model  for  this  case  is  now  exactly  the  same  as  for  the  previous  case,  a 
flat,  horizontal  bottom,  and  the  gravity  data  may  be  corrected  to  a  free-air 
(eq.  15)  or  Bouguer  (eq.  9)  anomaly  as  outlined  previously.   Grant  and  West 
(1965,  p.  240)  state  that  "unless  the  heights  (or  depths)  of  the  topographic 
features  exceed  one-twentieth  of  their  distance  from  the  reference  point,  they 
are  usually  omitted".   This  would  mean  for  this  case  that  the  terrain  correc- 
tion is  probably  unnecessary  for  a  bottom  slope  of  3°  or  less. 

To  correlate  the  gravity  values  taken  over  the  continental  shelf  from  the 
underwater  meter  with  sea-surface  gravity  measurements,  it  is  probably  best  to 
first  make  terrain  corrections  to  the  bottom  data  before  calculating  the  free- 
air  anomaly  using  eq.  15.   This  statement  is  based  on  the  concept  of  downward 
continuation  (Peters,  1949;  Trejo,  1954)  which  shows  that  for  a  relatively 
horizontal  bottom  the  terrain  effect  on  a  bottom  gravimeter  is  greater  than  on 
a  meter  at  the  sea  surface  at  some  distance  above  the  same  station.   In  shallow 
water  where  the  vertical  angle  from  a  given  portion  of  bathymetric  relief  to 
the  bottom  station  may  be  greater  than  the  similar  angle  from  the  feature  to 
the  surface  meter,  the  terrain  effect  may  be  greater  for  the  surface  meter, 
requiring  terrain  corrections  to  be  made  to  the  data  from  both  meters. 


Irregular  bottom 

Although  this  case  illustrates  the  usual  model  found  on  the  continental 
terrace,  the  previously  discussed  cases  are  applicable  approximations  in  many 
areas.   Along  the  central  California  coast  in  the  vicinity  of  Monterey  Bay  the 
presence  of  submarine  canyons  which  begin  just  seaward  of  the  shoreline  and 
the  rugged  coastal  mountain  range  cause  a  rather  narrow  continental  shelf  (about 
10  km  wide)  and  a  very  steep  continental  slope  (in  places,  bottom  slopes  of  up 
to  45°  have  been  found) .   Rock  outcrops  and  fault  scarps  on  the  shelf  cause 
many  bathymetric  irregularities.   The  model  for  this  case  may  be  broken  into 
the  same  three  segments  that  were  discussed  in  the  previous  case,  and  the  general 
method  for  making  terrain  corrections  is  the  same,  although  Sandberg's  correc- 
tion tables  for  an  inclined  plane  are  probably  not  appropriate.   Just  south  of 
Monterey  Bay,  Souto  (1973)  calculated  terrain  corrections  of  5.5  to  5.7  mgal 
for  stations  at  water  depths  of  60  to  75  m  near  the  Carmel  Submarine  Canyon, 
a  tributary  of  the  Monterey  Submarine  Canyon.   Inclinations  of  28°  were  found 
on  the  continental  slope  in  this  area.   With  Bouguer  anomalies  ranging  from 
+23  to  +30  mgal  for  this  area,  the  terrain  correction  is  an  important  factor 
in  data  reduction. 


CONCLUSIONS 

Correcting  underwater  gravity  data  presents  some  unique  problems  to  an 
investigator.   If  terrain  corrections  are  necessary,  it  appears  to  this  writer 
that  they  should  be  applied  first  to  the  observed  gravity  data  to  reduce  the 
bathymetry  to  a  model  of  a  flat,  horizontal  bottom.   Next  the  upward  attrac- 
tion of  the  water  above  the  station  should  be  removed  and  the  station  should 
be  relocated  on  the  reference  spheroid  with  a  free-air  correction.   Finally, 
water  should  be  again  added  to  the  ocean  beneath  the  spheroid  if  a  free-air 
anomaly  chart  is  desired  for  correlation  with  sea-surface  gravity  data,  or 
rock  should  be  added  if  a  Bouguer  anomaly  chart  is  required  for  correlation 
with  land  gravity  surveys.   This  order  of  application  of  corrections  should 
make  more  apparent  the  physical  meaning  of  each  step  and  aid  in  determining 
which  steps,  if  any,  may  be  simplified  by  approximation  for  special  bathymetric 
cases  or  eliminated. 
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